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Abstract 

We present a way of treating the problem of the interaction of a single trapped 
ion with laser beams based on successive aplications of unitary transforma- 
tions onto the Hamiltonian. This allows the diagonalization of the Hamil- 
tonian, by means of recursive relations, without performing the Lamb-Dicke 

approximation . 
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I. INTRODUCTION 



Trapped ions interacting with laser beams [ttlQ have become an extremelly interesting 
system for the investigation of fundamental physics, e.g. the generation of states of the har- 
monic oscillator [Q-|5|] as well as for potential applications such as precision spectroscopy 
and quantum computation J/J. The theoretical treatment of the interaction of a trapped ion 
with one or several laser beams constitutes a complicated problem, being the full Hamilto- 
nian highly nonlinear. Therefore approximations are normally required, such as, for instance, 
the Lamb-Dicke approximation, in which the ion is considered to be confined within a region 
much smaller than the laser wavelength. This makes possible to obtain Hamiltonians of the 
Jaynes-Cummings type, in which the center-of-mass of the trapped ion plays the role of the 
field mode in cavity QED. Recently ||, it has been suggested a new approach to this prob- 
lem, based on the application of a unitary transformation which linearizes the total ion-laser 
Hamiltonian. Moreover, in that approach || it is obtained a Jaynes-Cummings type Hamil- 
tonian including counter-rotating terms. A rotating wave approximation (RWA) may be 
performed provided we are in the special "resonant" regime, or Q « u, where Q is basically 
related to the laser intensity and u is the frequency of the ion in the magneto-optical trap. 
The RWA performed this way allows the exact diagonalization of the transformed Hamilto- 
nian, but on the other hand poses an upper limit to the possible values of the Lamb-Dicke 
parameter. In this paper we propose a method of diagonalization of the ion-laser Hamilto- 
nian entirely based on unitary transformations. As a first step we linearize the Hamiltonian, 
as it is done in |§ . The resulting Hamiltonian is further transformed, and cast in a form that 
it is suitable for diagonalization. Although the diagonalization procedure is an approximate 
one, based on recursion relations, it is not based in an expansion on any of the parameters 
of the problem. Therefore, the found solution does not depend on any specific regime, such 
as the Lamb-Dicke regime, for instance. 

This paper is organized as follows: In section 2 we obtain the linearized form of the 
full Hamiltonian. In section 3 we show how the resulting Hamiltonian may be further 
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transformed in order to be diagonalized. In section 4 we summarize our conclusions. 



II. LINEARIZATION OF THE HAMILTONIAN 



We consider a single trapped ion interacting with two laser (classical) plane waves (fre- 
quencies uji and U2), in a Raman-type configuration (u>l = uj\ — u 2 ). The lasers effectively 
drive the electric-dipole forbidden transition \g) <-» |e) (frequency ujq), with 5 = Uq — uj^. 
We end up with an effective two-level system for the internal degrees of freedom of the atom 
and the vibrational motion, after adiabatically eliminating the third level in the Raman 
configuration. This situation is described, in the atomic basis, by the following effective 
Hamiltonian [01 



H = h 



(1) 



being X = a + a\ h = a/ a, r] the Lamb-Dicke parameter, and a) (d) the ion's vibrational 
creation (annihilation) operator. 

By applying the unitary transformation 
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to the Hamiltonian in equation (HD, where D(/3) = exp(/3d^ — j3*a) is Glauber's displacement 
operator, with (3 = —ir]/2, we obtain the following transformed Hamiltonian 
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where Y = a — &t, and A = \t\v. This result holds for any value of the Lamb-Dicke 
parameter 7]. We have thus transformed our original Hamiltonian in equation ([[]) into a 
Jaynes-Cummings-type Hamiltonian, and therefore its exact diagonalization is allowed pro- 
vided we perform the rotating wave approximation (RWA). However this imposes limitations 
on the values of the Lamb-Dicke parameter. Because the "effective coupling constant" of the 



transformed Hamiltonian is given by A = \r\v, we shall have \r\v <C v if we want to neglect 
the counter rotating terms (RWA). Moreover, we still have to be "in resonance", or v ~ 2Q. 
In this regime the trapped ion-laser system suitably prepared exhibits long-time-scale re- 
vivals (superrevivals), as discussed in ||. Here we are not going to perform the RWA. We 
propose instead a way of further transforming the Hamiltonian in equation (|3|) in order to 
allow its diagonalization in an approximate, although nonperturbative way. In what follows 
we are going to present a novel method of transforming the Hamiltonian, showing how it 
may be diagonalized. 



III. DIAGONALIZATION OF THE HAMILTONIAN 



We have succeeded in transforming the ion-laser Hamiltonian in Equation (0) into a more 
tractable form (|3|). Nevertheless, its exact diagonalization has nor been achieved yet. The 
problem is either treated exactly after performing the RWA, or by means of approximate 
methods such as perturbative expansions or even numerically |TD[. Here we propose a 
different approach to that problem, based on unitary transformations. We are going to 
restrict ourselves to the resonant case 5 = 0. After discarding the constant term 



which just represents an overall phase factor, we obtain 

n = n 
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Now we define the following unitary transformations 
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These transformations have interesting effects when applied to the Jaynes-Cummings type 
Hamiltonian in (|j). First we apply Ti, which yields 



1 1 



Q AY 
AY -tt 



\( i A 

1 1 



V 



V 



-iXY n 
n iXY 



(7) 



i.e., the operators are transposed to the diagonal. Then we apply T 2 to the matrix which 
resulted form the operation above, such that 
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where it has been used the fact that (— l) n d(— 1)" = —a (the same for d^). We have now 
that Y is multiplied by the identity matrix. At this stage we apply T\ again, rearranging 
the terms in such a way that we obtain a Hamiltonian diagonal in the atomic state basis, or 
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Diagonalization of the total Hamiltonian becomes an easy task now, because the atomic 
part is already in a diagonal form. The total transformed Hamiltonian may be conveniently 
expressed in an operator form as 

n' = flfSflftHff&fi = h(yh + iXY - Qov(-l)") . (10) 

A general expression for the eigenstates of the Hamiltonian in expression flTU| ) is 

\*i) = \<fif)\9) + \ti)\e), (ii) 

with \iff) = Y] n Cni\n) and \(p 9 ) = X)„C^ z |n). From the eigenvalues equation H'l^fi) = 
Ai\^Ji) we obtain 



fa + iXY - n(-if) y: c^\n) = a,y: c:», 



(12) 



and 



(vh + zAY + U(-lf) £ C 9 Jn) =AiJ2 



(13) 
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The expansion coefficients {C's) may be obtained by means of recursion relations. For 
instance, for the coefficient Cf n . From equation ( [TT|) we may write 

E (nCS» +iAC^Vn|n - 1) - iAC^Vn + T|n + 1) - n(-l)»C^|n>) = A^C^n). 

n n 

(14) 

After rearranging some of the terms, we have the following relation between the coefficients 



By multiplying expression (|14|) by (0|, we also obtain the first term 

Ch = -i^Cl, (16) 

Similar relations may be found for the coefficients C 9 nl . Of course the normalization con- 
dition Yi, n |Cn«| 2 + |Cnzl 2 = 1 should be satisfied. Having diagonalized the transformed 
Hamiltonian, the evolution of the state vector becomes a trivial task. For that we have to 
express a generic state in terms of the basis states, or 

i 

The choice of a specific initial state determines the set of coefficients A\. For instance, if 
we initially prepare the trapped ion in the state |^(0)) = 4||/3) (\g) — |e)), the transformed 
state will read 

|V<0)) = ffflflF\*(0)) = |0)|e) = J2M*i)- (18) 

i 

The coefficients A\ may be determined by multiplying the left hand side of Equation flI8|) 
by (n\(e\ and (n\(g\, so that we obtain a set of coupled equations for the coefficients A[. It 
is now easy to calculate the time-evolution of the state vector, namely 

\ij(t)) =exp(-iH't/h)J2M*i) ='£A l exp(-iMt/h)\V l ). (19) 

i i 

The next step is to apply on \ip(t)) the sequence of transformations TiT^T/T^ backwards 
in order to recover the wanted solution for the state vector which describes the trapped ion 
system interacting with the laser beams. 
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IV. CONCLUSIONS 



We have presented a novel way of treating the problem of the interaction of trapped 
ions with laser beams entirely based on unitary transformations. The system Hamiltonian 
is successively modified and cast in a more tractable form. We propose a method of di- 
agonalization of the transformed Hamiltonian by means of the construction of recursive 
relations for expansion coefficients in the Fock state basis. Despite of not being an exact 
diagonalization, there is no need of performing any approximation, such as the Lamb-Dicke 
limit, for instance. Interesting possibilities for the investigation of such a system in different 
regimes are opened up, given that there are no restrictions on the relevant parameters and 
no approximations have been made so far. 
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